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^S^ • Abstract: The irreducible background from Z{-^ i/z^)+jets, to beyond the Standard 

II . Model searches at the LHC, can be calibrated using 7+jets data. The method utilises 

^^, the fact that at high vector boson pT (3> Mz), the event kinematics are the same for the 

,^ I two processes and the cross sections differ mainly due to the boson-quark couplings. The 

method relies on a precise prediction from theory of the Z/j cross section ratio at high pT, 

Cn ■ which should be insensitive to effects from full event simulation. We study the Z/7 ratio for 

^f^ ' final states involving 1, 2 and 3 hadronic jets, using both the leading-order parton shower 

^^ ■ Monte Carlo program PythiaS and a leading-order matrix element program Gambos. 

00 ' 

^vq . This enables us both to understand the underlying parton dynamics in both processes, 

r — [ , and to quantify the theoretical systematic uncertainties in the ratio predictions. Using a 

^^ [ typical set of experimental cuts, we estimate the net theoretical uncertainty in the ratio to 

be of order ±7%, when obtained from a Monte-Carlo program using multiparton matrix- 
elements for the hard process. Uncertainties associated with full event simulation are 
L^ ■ found to be small. The results indicate that an overall accuracy of the method, excluding 

;_! ' statistical errors, of order 10% should be possible. 
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1 Introduction 

At the CERN Large Hadron Collider (LHC), Z{-^ i^i/)+jets is an important Standard 
Model (SM) background to new physics processes that give rise to missing transverse 
energy + jets signals. In principle, the related process Z{-^ e"'"e~,/i"'"/i~)+jets provides 
a way to calibrate this background, although in practice the number of such events may 
be too small to do this with sufficient precision. It has therefore been proposed [1-5] to 
use a related calibration process, 7+jets production, which has a much higher rate.^ The 
key point is that at high transverse momentum, pT ^ Mz, 7 and Z production are very 
similar; indeed the only expected difference in rate comes from the different electroweak 
couplings of photons and Z bosons to quarks, which are of course very well determined. 
There is also no branching ratio suppression for photon production. The measured rate of 
7+jets production, coupled with theoretical knowledge of the ratio of Z and 7 SM cross 
sections, can therefore be used to accurately predict the Z{-^ i/i/)+jets background. 

Although the Z and 7 cross sections have a simple theoretical relationship at high 
vector boson pT, care is needed when estimating the theoretical ratio of Z and 7 SM cross 
sections. Especially for more than one jet, there are matrix element contributions to the 
cross section that may not be included in parton shower Monte Carlos such as Pythia 
[7, 8] or Herwig [9]. In addition, it is important to quantify the theoretical uncertainty 
on the ratio (from the choice of PDFs, QCD scales etc.), since this will propagate through 
to the overall uncertainty on the background estimate. A detailed theoretical study of the 



^The same technique has also been used by the CDF collaboration [6] to estimate a similar background 
in weak boson pair production in a hadronic final state at the Tevatron. 



V+ jets (y = 7, Z) cross sections and uncertainties is therefore required, to supplement the 
information from Monte Carlo event simulation. Finally, it is important to establish how 
well the theoretical precision survives under conditions closer to the experimental analysis, 
e.g. including effects from full event simulation, jet reconstruction, detector acceptance, 
experimental cuts etc. 

In this paper we report on such a study. Since our main interest is in estimating the 
missing Et distribution in events with multijets, we focus primarily on the inclusive vector- 
boson (Z or 7) transverse momentum (pt) distributions, particularly at high. pT S> Mz- We 
first analyse the ratio of the Z+1 jet and 7+I jet pT distributions from a general theoretical 
perspective using a program for up to and including 3 jets based on exact leading-order 
parton-level matrix elements. We then reproduce these results using PythiaS [7, 8] (at 
parton level) . This establishes that the ratio of the two process cross sections is theoretically 
robust, particularly at high px- We predict the value of the cross section ratio using 
'typical' experimental cuts, and estimate its theoretical uncertainty. We then consider the 
corresponding 2- and 3-jet cross sections, again comparing the exact leading-order matrix 
element results with those obtained from PythiaS. Finally, we use our results to assess 
the systematic uncertainties on the missing transverse energy + jets background obtained 
from the photon + jets cross section using this method. The ratios predicted by the two 
alternative approaches, based on "pure" matrix elements or parton showers, are illustrated 
as well as used to constrain the systematic uncertainties related to the commonly used, and 
also more optimal, scenario where matrix elements are matched with the parton shower 
used in the MC simulation. In the following sections, V refers to the vector bosons Z or 7. 

In a recent paper [10], a similar study was carried out for the V + 2 jets cross sections 
using two approaches: next-to-leading order in pQCD applied at parton level, and exact 
leading-order matrix elements interfaced with parton showers (ME+PS) as implemented in 
Sherpa [11]. Agreement between our results and those of [10] confirms that over much of 
the relevant phase space, and particularly for inclusive quantities, the effect of higher-order 
corrections on the leading-order Z, 7 cross section ratios is small. In our study, we consider 
also the 1- and 3-jet ratios and, more importantly, we extend the analysis to hadron level 
using PythiaS. This allows for an analysis similar to the ones within the LHC experiments. 
We also examine the dependence of the ratios on the parton distribution functions, since 
these receive different weightings in the Z and 7 cross sections and therefore do not exactly 
cancel in the ratio. Of course ultimately one would wish to evaluate all these cross sections 
consistently at next-to-leading order, using the methods described in [10] for the 2-jet 
case. 

2 y+jets production in leading— order perturbative QCD 

In the Standard Model, the coupling of photons and Z bosons to quarks q are, respectively, 

— tc 

-ieQg-f^" and j>^(v ag-f5), (2.1) 

2 sm Ow cos 0w 

where Qg, Vq and Og, are respectively the electric, vector and axial neutral weak couplings 
of the quarks, and 9w is the weak mixing angle. For hadron collider processes such as 



gg — )• V + ng OT qg ^ qV +{n— l)g, both of which contribute to ^ + n jets production, the 
matrix elements squared will contain factors of Q^ or {v'l + ai)/Asm^ 6w cos^ 9w for 7 or Z 
respectively. The only other difference in the matrix elements comes from the non-zero Z 
mass^, which will appear in the internal propagators and phase space integration. Sample 
Feynman diagrams for F + 1 jet production are shown in figures l(a,b), and for V + 2 jet 
production in figures l(c-f). 
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Figure 1. Sample Feynman diagrams for F + 1, 2 jets production, where V = ^,Z. 



For high priY) (^ Afz) production we would therefore expect the Z and 7 cross 
sections to be in the ratio 
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Substituting sin Oyy = 0.2315, we obtain R^ = 0.906 and Rd = 4.673. In practice, of 
course, the cross sections will receive contributions from all quark flavour types, and so 
R = a{Z)/a{'y) will be a weighted average of the Ru and Rd values, i.e. 
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^For the purposes of this discussion, we treat the Z as an on-shell stable particle. In practice, Z decay 
will also form part of the matrix elements. 



in an obvious notation, where {u) and (d) are the typical values of the n~type and d-type 
quark parton distribution functions (PDFs) in the cross section. Figure 2(a) shows i? as a 
function of the ratio {d)/{u). We would expect that where large x values are probed, for 
example at very high priy), the ratio would approach the i?„ value since d{x)/u{x) — )■ as 
x — )• 1, see figure 2(b). For moderate px values at the LHC, {x) ~ 0.1, which corresponds 
to {d)/{u) ~ 0.6 and therefore R ~ 1.4. 

The simple connection (2.3) between the vector boson cross section ratio and the 
initial state quark fiavour is, however, broken for njets ^ 2. Consider for example the 
sample Feynman diagrams of figure 1(e) and (f). For the former 'four-quark' diagrams, 
the vector boson can be emitted off any of the external quark legs and so the numerator and 
denominator of the ratio R depend on more complicated products of quark distributions. 
Because at high x uu scattering will be relatively more dominant than dd scattering, we 
would expect that the value of R for such processes would be closer to Ru than to Rd, 
compared to the 1-jet ratio. On the other hand, for the 55-scattering diagrams, figure 1(f), 
the ratio of the corresponding cross sections is (ignoring the Z mass) R = ^ Zg/ ^ 7g, 
where the sum is over the final state quark (antiquark) flavours, and the dependence on the 
initial state (gluon) distributions cancels. By way of illustration, with 5 massless flavours 
we obtain R = 1.933. As we shall see below, the four-quark contribution is more important 
at high pt than the gg contribution, and the net effect is to reduce R slightly compared to 
the 1-jet case. 
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Figure 2. (a) Dependence oi R = a{Z)/a{'y) at the coupling constant level on the ratio of average 
d and u parton distribution values, see eq. 2.3. (b) u and d quark MSTW2008LO PDFS [13] at 
Q2 = 104 GeV^. 



Figure 3(a) shows the Z + 1 jet and 7 + 1 jet cross sections^ as functions of the vector 



'Expressions for the matrix elements can be found, for example in Chapter 9 of [12]. 



boson transverse momentum at ^/s = 7 TeV and 14 TeV. Standard PDG values of the 
electroweak parameters are used, and the PDFs are the leading-order MSTW2008LO set 
[13] with renormahsation and factorisation scale choice ^r = fip = PTiV), and the Z is 
treated as an on-shell stable boson. The acceptance cuts are |2/(^, i)| < 2.5 and pt{V,J) > 
40 GeV, where y is the rapidity. Figure 3(b) shows the ratio of the Z and 7 distributions. 
We see the expected behaviour of a roughly constant ratio at large priV) ^ Mz lying 
between the Ru and Rd values defined above. Although above pr ~ Mz the ratio does 
exhibit a plateau region, at very large px we begin to see a slight decrease, as the high-x 
behaviour of the d/u PDF ratio drives the ratio down towards the Ru value. At 14 TeV, 
the empirical laige-pT value of i? ~ 1.4 is consistent with {d)/{u) ~ 0.6, see figure 2(a). 
This in turn is consistent with u and d PDFs probed in the x ~ 0.1 region, see figure 2(b). 
At the lower collider energy (7 TeV), higher x values are sampled for the same pT, and the 
Z/7 ratio decreases slightly, moving towards the Ru value. Note that the ratio curves in 
figure 3(b) can be reasonably well approximated by 
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with n w 1.2, illustrating the expected Z mass suppression relative to the photon distribu- 
tion for Pt < Mz- 
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Figure 3. (a) Z and 7 pT distributions at ^/s ~ 7 TcV and 14 TeV, with parameters and cuts as 
described in the text, (b) Ratio of the Z and 7 pt distributions. 



3 Parton level analysis 



In order to study the ratio of the Z and 7 distributions in a realistic experimental environ- 
ment, we need to use an event simulation Monte Carlo program. We use PythiaS for this 



purpose. At the same time, we want to understand the difference between the scattering 
amplitudes embedded in PythiaS and the amphtudes obtained using exact QCD matrix 
elements for multijet production. For the latter we use the program Gambos, an adapta- 
tion of the Giele et al. Vecbos program [14] for W, Z + n jets production, in which the 
weak boson is replaced by a photon. 

First, we compare the PythiaS and Gambos results at the parton level. This serves 
to check the consistency of the results from the two programs, when configured as similarly 
as possible, and provides a common middle step between the matrix-element (ME) V+ jets 
results at parton level produced by Gambos and the results for fully simulated events from 
PythiaS. 

The PythiaS results are obtained using the LO (2 — )■ 2) processes, qq — )• Vg and 
qg — )• Vq where V = -f or Z, corresponding to the Feynman diagram types shown in 
figures l(a,b). Events with > 2 jets are generated by parton showering off the initial and 
final state partons. This means that processes such as those shown in figures 1(c) to 1(f) are 
included, albeit with an approximation to the exact matrix elements, and we can therefore 
expect differences between Gambos and PythiaS results for the Z/7 ratios with > 2 jets. 

In order to produce results directly comparable with Gambos, the following settings 
are used as default in PythiaS: 

• PDFs: MSTW2008LO; 

• Strong: as{M'^) = 0.13939, with one loop running; 

• EM: a£;A/(M|) = 1/127.918, with one loop running; 

• Weak: sin'^{ew) = 0.2315; 

• Scales: Renormalisation and factorisation scales, ///j = fip = Pt(^); 

• Rapidity, transverse momentum and separation cuts on the final-state Z, 7 and jets 
as described in the previous section. 

In the following sub-sections we first compare results for the y + 1 jet distributions, then 
discuss the theoretical uncertainties on the corresponding Z/7 ratio, and finally compare 
the y + 2, 3 jet results from the two programs. 

3.1 V + 1 jet results 

The LO matrix elements used in the PythiaS processes are the same as those used in 
Gambos for the 1^ + 1 jet case, and a parton level comparison between the two programs 
should therefore give identical results. Note that the V and the jet correspond to the 
outgoing partons of the hard process in PythiaS, without any further simulation of the 
event. The differential Z and 7 cross sections and their ratio predicted by PythiaS are 
shown in figure 4 for pp collisions at 7 as well as 14 TeV. The results show the same 
characteristic features already seen in the Gambos predictions in figure 3, i.e. the Z cross 
section, excluding any branching ratios, is smaller than the photon cross section at small 
Pt due to the mass suppression, but is roughly proportional to, and slightly larger than. 



the 7 cross section at pT ^ Mz- The ratios obtained from PythiaS and Gambos are 
compared in figure 5, for 7 TeV and 14 TeV colhsion energies. Evidently there is good 
agreement between the two programs, as expected. As seen in the plot, the photon cross 
sections from the two programs agree perfectly, whereas a small difference between the Z 
cross sections, <5%, is visible. This difference is due to the way in which the PythiaS 
generator treats the Z boson as a resonance, in contrast to Gambos where the Z is treated 
as a real particle, and this was confirmed by producing a Gambos like process in PythiaS 
which reproduced the same results. Since the Z boson is generally treated as a resonance 
in MC programs, this diff'erence is not considered as a source of uncertainty. 
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Figure 4. (a) DifFerential cross section as a function of the vector boson pr for the process 
pp ^ V + 1 parton with V = j, Z from PythiaS, and (b) the ratio of these. 



3.2 Theoretical uncertainties 

Having established numerical agreement between the two programs, we can use PythiaS 
to investigate the theoretical uncertainty on the Z/j ratio at high px- There are a number 
of sources of these, which we address in turn. 

First, we consider the dependence on the PDFs used in the calculation. As argued 
in section 2, the Z/'j ratio at high pT is sensitive to the d/u parton ratio at large x. To 
study the possible variation in this ratio, we investigate the spread from using the different 
eigenvectors of the MSTW2008LO set and we compare these predictions with those from 
two (older) leading-order PDF sets, CTEQ5L [15] and GRV98 [16], shown in figure 6. The 
latter should yield a conservative estimate of the PDF dependence. The impact on the 
Z and 7 distributions and their ratio is shown in figure 7. Note that the CTEQ5L and 
GRV98 PDFs give respectively softer and harder Z and 7 pT distributions, which is an 
artefact of the underlying quark and gluon PDF behaviour, see figure 6(a), but that the 
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Figure 5. Effect on the boson pT and the ratio from using a matrix-element generator (Gambos) 
and PYTHIA8 at 7 and 14 TeV. 



effect largely cancels in the ratio. The residual small differences in the cross section ratio 
can be understood in terms of the corresponding small differences in the d/u ratio for the 
various sets, shown in figure 6(b). We ascribe a conservative ±4% PDF uncertainty to the 
R{Z/^) ratio at high priV)- 

In addition to varying the PDFs, we use different choices for the renormalisation and 
factorisation scales, fiR and fip, in order to mimic the effect of higher-order pQCD cor- 
rections not included in either the PythiaS or Gambos analyses. In particular, we use 
multiples of the default scales fiji = fip = PTiV), and two variants of this: the arithmetic 
and geometric means of the final-state transverse masses in the 2 — )• 2 hard process, fi'^^.-^ = 
(m^., + ml^)/2 = {2pT{Vf + m{Vf)/2 and /i^^^ = mTimT2 = Pt(^)Vpt(^)2 + m(y)2. 
Note that for the photon, these scales are identical to the default scale pt{i)- Figure 8 
presents the corresponding impact on the differential cross sections, da/dpT, as well as the 
cross section ratio, R{Z/^). The results show that although the variations have significant 
effects on the differential cross sections, as expected, the Z/^ ratio remains stable in the 
regime pT ^ Mz, and for priY) > 100 GeV all variations of the ratio are within ±3%. 

Indeed the only sizable effect on the ratio related to the scales is observed from the 
different choices of the scale //r, which becomes visible at pT 'C Mz- Any choice of scale 



of the form upT will of course cancel in the Z/7 ratio, but scales of the form Ka/p"^ + My 
will give different results for low pT ~ My- The size of this effect was also shown to 
be consistent with the ratio as{nR{Z))/as{fiR{'y)) using the same one-loop formula as 
in PythiaS and Gambos. No similar effects are observed from different choices of the 
factorisation scale np, since the PDFs vary only weakly with the factorisation scale at the 
X values probed by these cross sections. 

Note that in the above analysis we have used the same form of scale variation simul- 



1.4 
1.2 
1.0 
O 0.8 
^ 0.6 
0.4 
0.2 
0.0 



■ 1 1 




- MSTW08LO 


- 


■ GRV98LO 


■ 


L CTEQ5L 


■ 


: Q' = 10'GeV' 


- 


l^^u^ 


-_ 


. . . 


\ ^ 





\^] 



0.01 



0.1 




(a) 



(b) 



Figure 6. (a) Comparison of up and down quark distributions from MSTW2008LO, CTEQ5L and 
GRV98. (b) d/u ratios compared to MSTW2008LO. 



taneously in both the numerator (Z) and denominator (7) cross sections. As pointed out 
in Ref. [10], this gives a much smaller scale variation than if the scales are varied inde- 
pendently in the two cross sections. However, we argue that if we select Z and 7 events 
for which the kinematics of the (colour-singlet) vector bosons and the jets are the same, 
and if the energies and momenta are large enough such that the Z mass can be neglected 
(e.g. pt ^ Mz), then the higher-order pQCD corrections to both cross sections should 
essentially be the same and should therefore largely cancel in the ratio. 

3.3 y + 2, 3 jets results 

For the V + 2 jets production cross sections there is an additional complication in that 
the high-pT^ photon can be emitted collinearly to a high-p^ quark, with the transverse 
momentum of the pair being balanced by an 'away side' quark or gluon. For massless 
photons, quarks and gluons the matrix element is singular in this configuration and so the 
closer the photon is allowed to approach the quark, the smaller the Z/7 ratio becomes, 
since there is no such collinear singularity for Z production. To regulate the singularity, we 
impose a A.R{V, j) > AR^i^ isolation cut, where V = ■jjZ and j = q, g.^ We also, of course, 
need to impose rapidity, transverse momentum and jet-jet separation cuts on the quark 



*Note that the requirement of photon isolation becomes more subtle beyond leading order in pQCD, 
since in this case the partonic jets can have non-zero 'width' and the choice of jet algorithm influences the 
analysis. This issue is addressed in detail in Ref. [10]. In our case we will be studying photon isolation for 
the full PythiaS event simulation including hadronisation and experimental cuts. Note also that we neglect 
contributions to the 7 cross section involving photon fragmentation functions, i.e. /'~*^(2, Q^). With the 
strong isolation requirements and high transverse momentum values used in our study, we expect such 
contributions to be small. 
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and gluons jets. For illustration, we choose ptO) > 40 GeV, |y(j)| < 2.5 (as for the 1-jet 
study above), and AR{j,j) > 0.4 to represent 'typical' experimental cuts. Figure 9(a) 
shows the ratio of y + 2 jet cross sections, as calculated using Gambos, for different 
values Ai?min = 0.05,0.1,0.2,0.4,0.6 at 7 TeV. The ratio at high pT shows the expected 
dependence on the minimum separation. Note also that the ratio becomes insensitive to the 
isolation cut when the minimum separation becomes large, since far from the singularity 
the Z and 7 phase space are affected more or less equally. A similar dependence on AiJmin 
is observed for the 1^ + 3 jet ratios. From now on we take AR^i^ = 0.4 as our default choice 
and attribute a ±5% uncertainty of these results based on the difference with respect to 
AR^in = 0.6. 

The breakdown of the Z + 2 jet Gambos cross section at 7 TeV into the different 
subprocess contributions is shown in figure 9(b). We define these to be qq — )• Vgg scattering 
(e.g. figure 1(c)), qg,gq -^ Vqg (e.g. figure 1(d)), gg ^ Vqq (e.g. figure 1(f)) and qq -^ Vqq 
(e.g. figure 1(e)), where a sum over quarks and antiquarks is implied. Note that quark- 
gluon scattering is by far the most dominant in the kinematic region studied here, and 
its fractional contribution is roughly independent of priV)- The results also show that 
the second largest contribution in the 2-jets case comes from the qq subprocess, which 
approximately amounts to 20%. This is in contrast to the l~jet case, which is dominated 
by qg and qq scattering. The corresponding subprocess breakdown of the 7 + 2 jet cross 
section is similar. 

In figure 10(a) we show the Z/7 + 1,2,3 jet Gambos cross section ratios as a function 
oi pt{V), with Ai?inin = 0.4 and other cuts as before. We see that the 2,3 jet ratios are 
slightly smaller than the 1 jet ratio at moderate and high px- The small difference arises 
from three effects: (i) the dependence of the 2,3 jet cross sections on Ai^min = 0.4, (ii) the 
additional qq and gg scattering diagrams, the net effect of which is to decrease the 2,3 jet 
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Figure 8. Effects on (^,7) differential cross sections and cross section ratio after varying the 
scales {fJ.R-, (J-f)- The different fi scales are defined in the text and the denominators correspond to 
using the default scale choice, i.e. = {pt,Pt)- 



ratio, as already explained in section 2, and (iii) the fact that for a fixed pT, increasing 
the number of jets increases the overall invariant mass of the final-state system, and also 
therefore the values of the parton momentum fractions. This in turn decreases the d/u 
ratio, and also the Z/j ratio, see figure 2(a). 

For more than one jet, the additional jets in the PythiaS simulation are produced by 
parton showering. In order to characterise the difference with respect to exact MEs, the 
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Figure 9. (a) Ratio of the Z + 2 jets and 7 + 2 jets pt(I^) distributions at 7 TcV, for different 
values of the Ai?(y, j) > Ai?niin isolation cut. (b) Breakdown of the Z + 2 jet Gambos cross section 
into the different subprocess contributions defined in the text. 
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Figure 10. (a) Z/j + 1,2,3 jet Gambos cross section ratios as a function oi pxiV), with the 
isolation cut Ai?min = 0.4 . (b) Effect from multijet requirement on the PythiaS ratio (upper) 
and difference with respect to Gambos (lower). 



following events were used for the 2- and 3-jet results from PythiaS at the parton level, 
where the aim is to produce events with a similar parton level topology to the Gambos 
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events based on the multijet MEs. In the 2-jet case, events were generated with both 
initial (ISR) and final state radiation (FSR) enabled and only those with exactly 2 jets 
within the allowed pT and y acceptance were considered. These jets either correspond to 
the parton from the hard scatter together with an ISR emission or from events with no 
accepted ISR emission, but where the hard parton, due to FSR, branches into two partons 
within the allowed acceptance. In addition, only events with the required V— jet and 
jet-jet AR separation were considered. The same event selection was used for the parton 
level 3-jet case, where the three jets either correspond to the hard parton together with 
two ISR emissions, or a FSR branch of the hard parton together with one ISR emission, 
or two FSR emissions. The multijet ratios from PythiaS were found to be very similar 
to the one jet case. This is illustrated in the upper plot of figure 10(b), which shows the 
PythiaS -R2jet/-Rijet and Rsjet/Rijct ratios.^ The corresponding Gambos Z/7 ratios are 
therefore slightly smaller as shown in the lower plot of figure 10(b). These results illustrate 
the difference between the two individual approaches in the multijet case as well as the 
importance of using a multiparton ME based program in the actual analysis, where the 
precision related to the ME calculation was discussed above in connection with figure 9. By 
comparing the Gambos and PythiaS ratios at high pj-, we can extract correction factors 
for the PythiaS ratios to take account of the missing contributions, however, as seen in 
figure 10(b) these correction factors are not large. 

4 Full event simulation 

The ability of PythiaS to simulate full events was used both to investigate the robustness of 
the Z/'j cross section ratio as well as its potential use in estimating the Z ^ vv background 
in searches for new physics at the LHC. For simplicity, the experimental aspects related 
to the photon analysis attempt to follow as closely as possible what is commonly used 
in ATLAS analyses [17, IS], and for the new physics scenario we focus on the SUSY 
zero-lepton search [1, 2], where SM Z — t- z^P production is one of the main backgrounds. 
Due to the phenomenological nature of this study, we neglect any experimental photon 
inefficiencies, apart from the isolation criteria discussed below, as well as any backgrounds 
(e.g. TT*^ — )• 27), which in any case are expected to be relatively small in the high pT region 
of interest. 

4.1 Effects on the ratio 

The same PythiaS processes were used as in the parton-level study, but with the full 
parton shower, hadronisation, multiple interaction and particle decay simulation enabled. 
The default settings of vS.150 were used, for which general performance results can be found 
in [19]. The same selection was used as for the 1-jet parton-level results, here corresponding 



^The large deviation in the first bins is an artefact of the kinematic acceptance. In the 1-jet case, events 
will populate these bins according to the normal underlying distribution. However, since shower emissions 
can only occur with pT smaller than the hard process, the 2,3 jet events will be peaked toward the upper 
bin edge. Due to the sharp rise of the ratio at low pt, the average ratio value in these bins will therefore 
be significantly different for the higher jet multiplicities. 
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to an inclusive jet selection. The main differences with respect to the parton-level results 
come from using the final state boson momentum as well as from using jets reconstructed 
from the final state particles, rather than being represented by single partons. The jets were 
reconstructed using the FastJet library [20-22] and were based on all final-state particles 
except leptons and any photons with pj, > 30 GeV. The anti-fcj algorithm was used with 
a R parameter of 0.4, also in accordance with the ATLAS analysis. In the following, these 
results from full event simulation are referred to as obtained at particle level. 
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Figure 11. (a) Transverse energy inside the photon isolation cone as a function of the photon pr- 
(b) Z/7 cross section ratio as a function of the photon pT- Parton (part) and particle level results 
(reco) are shown as well as for events passing the isolation criteria (reco + iso). 



The experimental analysis in ATLAS uses a photon isolation criterion in order to sup- 
press QCD background and this quantity was found to be well described at MC generator 
level. This isolation criterion requires the transverse energy within a Ai?[= yAj/^ + Ac^] = 
0.4 cone around the photon {E^°) not to exceed 4 GeV. Since the anticipated use of the 
ratio here is to estimate the number of Z events from measured 7 events, the effect from 
this isolation requirement is also addressed. The fact that only isolated photons are con- 
sidered implies an even stronger photon-jet separation than the one used above for the 
parton level results, ensuring an acceptance where Gambos provides robust calculations. 
The mean transverse energy within the photon isolation cone is shown in figure 11(a) as 
a function of the photon pT- This plot is based on events which pass the above selection 
applied to the final state photon as well as the reconstructed jets. A small increase with pT 
is shown, but with values well below 4 GeV over the whole range. In spite of the increasing 
hadronic recoil with larger boson pT, a relatively constant inefficiency of about 5% was 
found over the full px range. 

Figure 11(b) shows the Z/7 cross section ratio at parton (part) and particle level (reco) 



14 



as well as at particle level where the isolation requirement is applied (reco+iso). Good 
agreement is evident between the results obtained at parton and particle level, where the 
difference is well below 5% at high p^", and the increase of the ratio due to the photon 
isolation criterion is of order 6% at high px- 

4.2 Background estimate for a zero lepton SUSY search 
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Figure 12. (a) Boson pT and (b) 77 distributions from events passing the 2-jet SUSY selection. 
Results from photon events both with (sel) and without (no sel) applying the photon analysis 
selection are shown together with results from Z events (Z). 



This section demonstrates the estimation of Z — t- i/p background for a zero lepton 
SUSY search using photon events. This is done using the PythiaS results and is meant 
to serve as a general example, since the method could be used also for other new physics 
searches where Z ^ vv production contributes with a significant background. The method 
involves the following steps: 

• Photon event selection. Select a photon event sample using a loose enough selection, 
with respect to the photon and jets, to contain as many events as possible that will 
pass the final selection. This is represented here by the criteria pri'y) > 45 GeV, 
|?7(7)| < 2.37, excluding 1.37 < |?7(7)| < 1.52 and E'^" < 4 GeV, based on the 
ATLAS photon analysis [17, 18]. 



• SUSY event selection. Apply the SUSY selection to the photon events, where the 
photon pt represents the missing transverse energy from the Z in the events to be 
estimated. A 2-jet as well as 3-jet SUSY selection is used, based on the ATLAS 
search. 2-jet (3-jet): priji) > 120 GeV, pt{J2) > 40 GeV, {prUs) > 40 GeV), 
|r?(j,)| < 2.5, pt{V) > 100 GeV, A(^(U,ii) > 0.4, PT{V)/m^jf > 0.3 and m^jj > 500 
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Figure 13. (a) Boson pT and (b) rj distributions from events passing the 3~jet SUSY selection. 
Results from photon events both with (sel) and without (no sel) applying the photon analysis 
selection are shown together with results from Z events (Z). 



GeV. Here ji represents the i leading jet and meff is the SUSY discriminating 
variable used in [1, 2], defined as the scalar sum of the pT from the jets and the 
boson in the event. 

• Subtract backgrounds and correct for experimental efficiencies. This is represented 
here only by the isolation efficiency. 

• Convert photon events, inside the acceptance of the analysis, to Z,yy events using the 
cross section ratio, R{pj-) ■ Br{Z — )• vv). As discussed in the previous sections, in 
an analysis of real LHC data, R{pj-) should be based on results from exact multijet 
MEs and using an appropriate jet selection. 

• Correct for acceptance constraints implied by the photon analysis, e.g. the 77(7) se- 
lection criteria. 

The main intention with this method is that all necessary corrections as well as theoretical 
input are related to the vector bosons, whereas all requirements with respect to the ex- 
perimentally more challenging reconstructed jets, are identical for the Z and 7 events. As 
shown in the previous section, the ratio is affected by requiring jets. However, since this 
effect is a consequence of changing the mixture of couplings imposed by the relevant initial 
partons together with their PDFs, it is small even for drastically different jet criteria and 
should be yet smaller with respect to experimental jet uncertainties, such as energy scale 
and resolution. 
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The photon event selection imposes some unavoidable criteria which are not experi- 
enced by the Z^y background and therefore has to be corrected for. The implications of 
this selection on the final sample are, however, relatively mild for the following reasons. 
The photon p-p requirement is significantly softer than the subsequent selection. A photon 
isolation criterion, of some kind, is required in order to obtain accurate calculations of 
the ratio. However, this is often also used in the SUSY selection for more experimental 
reasons, such as to prevent fake missing transverse energy caused by a high p-p jet, i.e. the 
A(/)(y, jj) requirement above. In addition, the fact that the Z and 7 processes have the 
same phase space when pj- ^ Mz means that at high p^^, the r/ distributions from the 7 
and Z events converge toward the same distribution, which becomes increasingly central 
with higher boson "px- Therefore any acceptance corrections will become the same for the 
Z as for the 7 events. 

The SUSY selection, based on the "PT of the boson and jets in the event, is then identical 
for the two event types and should not require any related corrections. Due to the fact 
that the bosons are recoiling against the hadrons, the SUSY selection will in principle act 
as a non-trivial high boson pr criteria. For this reason the 7 events passing the SUSY 
selection can be converted into Z events based only on the boson kinematics. Again due 
to the convergence of the Z and 7 phase space at high p^^, the cross section ratio becomes 
insensitive to the particular r\iy^ criteria used and is hence determined by the Pt(^)- The 
precision of this method is therefore mainly related to the photon analysis part, which is 
expected to be precise at high photon pTj and the theoretical knowledge of the Z/7 cross 
section ratio. 

In figures 12 and 13 the boson •p'x and 77 distributions are shown after the 2-jet and 
3-jet SUSY selections respectively. The distributions for photon events, with (sel) and 
without (no sel) passing the photon selection, as well as Z events (Z) are shown. The pT 
distributions show that both the SUSY selections mainly select events with a boson pT 
in the range 250 to 300 GeV. The ?/ distributions also show the similarity in shape of the 
distributions from the Z and 7 processes. Both the pr and r/ distributions show that the 
effect from the photon selection, i.e. acceptance and isolation requirements, is relatively 
small also after the SUSY selection and the difference between the Z and 7 results is 
consistent with the cross section ratio at relevant boson px values. In the r\ distributions, 
the difference from applying the photon selection reflects the impact of the isolation criteria 
alone and it is shown that the additional jet requirements from the SUSY selections do not 
change the isolation efficiency dramatically. The results from the 2-jet and 3-jet selections 
do have slight differences, but the overall characteristics discussed are the same. Figure 14 
shows the Z pT distributions from 7 events passing the two SUSY selections, converted into 
Z ^ vv according to the method outlined above, together with the results obtained from 
simulating Z ^ vv directly. Since the PythiaS ratios do not show any jet multiplicity 
dependence, the ratio shown in figure 11 (reco) was used in this analysis example for both 
the 2-jet and 3-jet results. As expected, the two distributions agree within the statistical 
uncertainty of the simulation, which is smaller than 5% in the bulk of the distributions. 

The uncertainties on the final background estimate, related to the experimental aspects 
of this analysis, are not covered by this study. However, as discussed above, the design 
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Figure 14. Differential cross section as a function oi pt{Z) for Z ^ vv events passing the (a) 
2-jet and (b) 3-jet SUSY selection. The predictions using 7 events (Zi,y from 7) are compared to 
results from direct MC simulation of the Z ^ vv process iZ^t,). 



of the method should limit the exposure mainly to the corrections from the photon event 
selection, which are expected to be precise at high boson px- In addition, as discussed in 
section 4.1, the impact associated with full event simulation on the cross section ratio was 
found to be small. As shown in figure 11(b), the ratio increases, by about 6% in the high 
•pT region, when including the isolation efhciency, whereas a significantly smaller effect is 
seen when going from parton to particle level. The final uncertainties on R[Z/^) from 
such effects are therefore expected to be at the percent level. A theoretical uncertainty, 
here with respect to the hard process calculations behind R{Z/^) ■ Br{Z — t- i^i^), that 
significantly exceeds 10% is therefore likely to dominate the uncertainty of the method. 

Table 1 presents the overall R ■ Br values obtained for the results shown in figure 14. 
This includes both the value from PythiaS {R ), that was actually used in the plot, as 
well as the corrected value [R ) based on the ME results. The table also includes the 
uncertainties associated with the ME calculation (eme), the scales (e^) and the PDFs 
(^pdf), see section 3. The results show a total uncertainty of ±7% and indicate that the 
theoretical uncertainty for results obtained by an appropriately configured MC program, 
which uses ME amplitudes for the hard jets, should be within 10%. The SUSY selection 
used in the experimental analysis will evolve with an increasing amount of LHC data. 
However, such an evolution is expected to effectively imply a harder boson pT requirement 
and since the uncertainties in table 1 are valid for pt(^) > 100 GeV, the same conclusions 
should hold also for harder selections.^ 



'At least up to pt{V) — 800 GeV, which is the maximum value included in this study. 
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Selection 


RP-Br 


j^ME . ^^ 


£me 


e^ 


£PDF 


£Tot 


2-jet 


0.254 


0.234 


5% 


3% 


4% 


7% 


3-jet 


0.246 


0.207 


5% 


3% 


4% 


7% 



Table 1. The overall i? • i?r values obtained after the SUSY selections, directly from PythiaS {R^) 
as well as corrected values with respect to the ME results (i?*^-^). The uncertainties associated 
with the ME calculation (emb), the scales (e^) and the PDFs {epdf) are also shown. 



5 Summary and conclusions 

One of the best methods to calibrate the irreducible background from 2'(— )• z/;y)+jets, to 
beyond the SM searches at the LHC, comes from using 7+jets data. The method utilises 
the fact that at high boson pT (^ Mz) the event kinematics converge for the two processes 
and the cross sections differ mainly due to the boson couplings. The advantage comes from 
large statistics, compared to alternative methods using Z{^ ee,/i/i)+jets events, together 
with the clean signature, with respect to experimental efficiencies and background, at high 
photon pt- Hence, a precise prediction from theory of the Z/j cross section ratio, R{Z/^), 
is required. The similarity between the two processes should allow for a robust prediction 
of R{Z/j), given careful attention to the modelling of the jets. 

The general dependence of R{Z/^) on the mixture of boson couplings, which is de- 
termined by the initial state partons of the relevant amplitudes and their corresponding 
PDFs, has been illustrated. Relatively accurate values can be obtained even using rough 
approximations in the 1-jet case, whereas a larger set of amplitudes becomes necessary 
when 2 or more jets are required. The ratios have been studied at parton level using both 
the (LO) PythiaS as well as the (multijet ME) Gambos programs, which allows us to 
disentangle effects associated with the two approaches. The impact from exact MEs when 
requiring different numbers of jets was found to be significant, but uncertainties were found 
to be within 5% for the acceptance from typical experimental cuts. The corresponding un- 
certainties related to the PDFs and scale choice were found to be less than 4% and 3% 
respectively. 

The PythiaS MC program was used to investigate effects on R{Z/'y) associated with 
full event simulation as well as performing a proof of principle analysis example. The 
effects investigated were found to be small and indicate that a theoretical precision^ at the 
10% level is required, in order not to significantly degrade the performance of the method. 
The total theoretical uncertainty was found to be 7%, indicating that the results obtained 
by MC simulations, including exact multijet MEs, should be within the 10% level. These 
results should also hold for similar 2- and 3-jet selections, given that the effective pT 
requirement on the boson is harder than what is used in the example analysis. Note that 
all our theoretical cross sections are evaluated in leading-order pQCD. It will be important 
to check, using for example the techniques of [10], that the ratio predictions are indeed 
stable - at least to the required accuracy - with respect to higher-order pQCD corrections. 



^Again, referring to the hard QCD process calculations behind R ■ Br. 
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Finally, one type of correction that has not been included in our theoretical study is 
high-order electroweak corrections. Although these are intrinsically small, and many will 
again cancel in the Zj^ ratio, there is an important class of correction involving W and 
Z virtual exchanges that does not cancel in the ratio. The impact on both the Z and 7 
distributions have been studied in [23, 24]. It was shown that non-cancelling Sudakov-type 
logarithms ~ a log (pt(^)^/?^vv') appear at highpr(y), and decrease the Z/7-I-I jet ratio 
by 6% (11%) at pt(^) = 300 (800) GeV [24]. However care is needed in the interpretation 
of this result, since the emission of real W bosons is expected to compensate the virtual 
Sudakov logarithms to some extent [25, 26]. It is therefore important to carry out a full 
analysis of higher-order electroweak corrections for the multijet processes and acceptance 
cuts studied in this paper. 
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